Stationary entanglement achievable by environment induced chain links 
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We investigate the possibility of chaining qubits by letting pairs of nearest neighbours qubits dissipating into 
common environments. We then study entanglement dynamics within the chain and show that steady state 
entanglement can be achieved. 
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I. INTRODUCTION 

It is nowadays well established that entanglement repre- 
sents a fundamental resource for quantum information tasks 
Ql]. However, being a purely quantum feature it is fragile 
with respect to enviromental contamination. Notwithstand- 
ing that, the possibility to achieve entangled states as station- 
ary ones in open quantum systems has been put forward in 
many different contexts (for what concern qubits systems see 
e.g. Refs.)|2, S]). The subject has attracted a lot of attention 
up until a recent striking experiment on long living entangle- 
ment ||4(] . The works on this topic can be considered as falling 
into two main categories: one where all qubits are plunged in 
the same environment |2] and the other where each qubit is 
plunged in its own environment 0]. In the former case the 
environment can provide an indirect interaction between oth- 
erwise decoupled qubits and thus a means to entangle them. 
In the latter case a direct interaction between qubits is needed 
to create entanglement, and usually to maintain it one has to 
also exploit other mechanisms (state resetting, driving, etc.). 

Here we consider a hybrid situation as depicted in Fig.l. It 
represents a sort of spin-^ chain dimerized by environments. 
In practice each environment induces a chain link between 
contiguous qubits. Hence, we can expect that a simple dis- 
sipative dynamics in such a configuration is able to establish 
entanglement along the chain without the need to exploit any 
other mechanism. Actually, we will show for the case of three 
qubits the possibility of achieving stationary entanglement for 
each qubits pair. The amount of entanglement results strongly 
dependent on the initial (separable) state. Also the depen- 
dance from the chain boundary conditions (open or closed) 
will be analyzed as well as a left-right asymmetry in qubit- 
environment interaction. 

The layout of the paper is the following: in Section II we 
introduce the model relying on physical motivations and we 
discuss the general dynamical properties; in Section III we 
restrict our attention to the three qubits case and investigate 
the entanglement dynamics in the open boundary condition; 
in Section IV we analyze the same system but with closed 
boundary conditions. Concluding remarks are presented in 




FIG. 1: Environments (ellipses) inducing chain links between con- 
tiguous qubits (gray circles). 



Section V. 



H. THE MODEL 

The model of Fig. 1 can be motivated by physically consid- 
ering two-level atoms inside cavities connected by fibers 
In such a scheme each atom-qubit can be thought as exchang- 
ing energy with the optical modes supported by the fiber. In 
turn this latter can modeled as an environment. Thus each 
qubit dissipates energy through two environments (one on the 
left and the other on the right). It happens that two contiguous 
qubits dissipates energy into the same environment. Then this 
environment mediates the interaction between the contiguous 
qubits. 

More specifically, let us consider at the ?'th site of a chain 
a qubit described by ladder operators af satisfying the usual 
spin-i algebra [af ,&~] = of. Let us also consider at the 
ith site of a chain radiation modes described by ladder opera- 
tors bij, b\ j satisfying the usual bosonic algebra [bij, b\ ■,] = 
Sjjr . Then, the interaction Hamiltonian reads 

^ = E E <V (''< i... + ki) + h-c (1) 
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By considering the b\ t j as environment's operators for the ith 
qubit, we can use standard techniques |0] to arrive at the fol- 
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lowing master equation 

dp 



(2) 



where {, } denotes the anti-commutator and we have assumed 
unit decay rate. 

Since we are interested on the steady state we have to no- 
tice that, given a master equation written in the standard Lin- 
bladian form, 



dt ^ 



LipL] - 
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(3) 



the uniqueness of the stationary solution is guaranteed if the 
only operators commuting with every Lindblad operator L\ 
are multiples of identity QV|] . 

In the case of Eq.© the cK~s commute with Lindblad oper- 
ators. Hence the steady state may not be unique, that is it may 
depend on the initial conditions. Due to that we need to study 
the full dynamics of the system. 



III. THE THREE QUBITS CASE WITH OPEN BOUNDARY 
CONDITIONS 

We restrict our attention to a chain of three sites. We first 
consider open boundary conditions. Then, the dynamics will 
be described by a master equation that can be easily derived 
from Eq.© 

^ = 7 [2(a 1 -+a 2 -)p(a+ + a+) 

-{(at + <t+ ) +a 2 - ),p}] 

+ (l-7)[2(^~+a 3 -) i 5( ( 7 2 + + ( T+) 

-{(^+°i){°2 +&s),f>}] ■ (4) 

Here we have considered the possibility for each qubit of hav- 
ing an asymmetric decay rate on the left and right environ- 
ments. This has been accounted for by the real factors 7 and 
(1— 7) with the assumption < 7 < 1. Clearly the symmetric 
situation is recovered when 7 = 1/2. 

By arranging the density matrix (expressed in the com- 
putational basis {|e), |<?)}® 3 ) as a vector v (e.g. writing 
Pi.j = v 8(i-i)+jX me master equation (0]i can be rewritten 
as a linear set of differential equations 

v(t) = Mv(t), (5) 

where M is a 64 x 64 matrix of constant coefficients given by 

M = ^(2L 1 ®L 1 -L[L 1 ®h-h®L\L 1 ) 

+ (1 - 7)(2i 2 ® L 2 - l\l 2 ® I s - I 8 ® l\l 2 )16) 

where 

(7) 
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with /„ the n x n dimensional identity matrix and a + = 

0_) ' CT = ^ 1 J ' ^^ en ' t ' le set °^ differential equa- 
tions <(5j> can be converted into a set of algebraic equations via 
the Laplace transform, v(s) = J °° exp(— st)v(t), i.e. 



sv(s) - v(0) = Mv(s). 



(8) 



Decoupling these equations one finds that the Laplace trans- 
forms of the density matrix elements are rational functions 
of polynomials and the inverse Laplace transformation can be 
performed analytically. The results are not explicitly reported 
because the expressions are too much cumbersome. 

Having the density matrix of the system, we can study the 
entanglement dynamics for each qubit pair of the system. To 
quantify the amount of entanglement between each of the 
qubits we use the concurrence HI]. We recall that to find the 
concurrence of a bipartite system described by the density ma- 
trix p, the following steps should be done: 

1) Find the complex conjugate of the density matrix in the 
computational basis and denote it by p*. 



2) Define p 



(<Jy®<Jy)p* (oy®a y ), where a y — i(a 



3) Find the square root of the eigenvalues of pp and sort 
them in decreasing order: {Ai, A 2 , A3, A4}. 

4) The concurrence is given by 

C = max{0, Ai - A 2 - A 3 - A 4 }. (9) 

A. Entanglement dynamics 

Figure|2]shows the evolution of entanglement between each 
qubit pair for |eee) initial state. As it can be seen in this fig- 
ure, when all the qubits are initially in excited state it takes 
longer time for the first and the third qubits to become entan- 
gled compared to the time needed to generate entanglement 
between the first and second or second and third qubits. As 
a consequence for not nearest neighborhood qubits we have a 
sudden birth of entanglement, i.e. it suddenly becomes non 
zero at times greater than zero (it does not smoothly increases 
starting from initial time) Q9D . 

If we start with different initial states, the entanglement be- 
haves differently in time. Figure [3jleft shows the value of 
concurrence in time for \ eeg) as initial state. Like the previous 
case, it shows that the time needed by the first and third qubits 
to become entangled is longer than the time needed by nearest 
neighborhood qubits. Moreover, we still have the entangle- 
ment sudden birth phenomenon, but this time it manifests not 
only for distant qubits (first and third) but also for the first and 
the second qubits. The interesting point in this figure is that 
the entanglement generation between the nearest neighbors is 
quicker if they are initially prepared in \eg) rather than |ee) 
state. The other possibility with two number of excitations 
in the initial state is \ege), for which the time evolution of 
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FIG. 2: Qubit pair concurrence vs time for \eee) initial state. Dashed 
line: C 1 2 , solid line: C 2 3 and dotted line: C 1 3 . 
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FIG. 4: Qubit pair concurrence vs time for \egg) (top) and 
(bottom) initial states. Dashed line: C 1 2 , solid line: C 2 3 and dotted 
line C 1 3 . 
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FIG. 3: Qubit pair concurrence vs time for \eeg) (top) and \ege) 
(bottom) initial states. Dashed line: C 1 2 , solid line: C 2 3 and dotted 
line: C 13 . 



generated between each qubits pair from the beginning, no 
matter how far they are from each other, i.e. we no more have 
entanglement sudden birth phenomenon. 

Finally, in the case of ggg) initial state there is no entangle- 
ment at any time because this state represent a fixed point of 
the Liouvillian superoperator at right hand side of the master 
equation (|4]i, or in other words Mv = with Vj = <5j,64- 



B. Stationary entanglement 

Taking the limit t — > oo in the density matrix elements, we 
arrive at the following general form for the steady state 



entanglement in shown in Figure[3]-right. This time the entan- 
glement sudden birth phenomenon only manifests for distant 
qubits (first and third). 

A big difference appears if the number of excitation of the 
initial state, reduces to one. Figures [4}left and-right show en- 
tanglement evolution for the initial states \egg) and \geg), re- 
spectively. As it can be seen in this figure, entanglement is 
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FIG. 5: Steady state entanglement versus 7 for different initial states. 

where / e K is determined by the initial state of the system. 
In particular we have the following correspondence: 
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To find the amount of entanglement in each qubits pair at the 
steady state, we first write the reduced density matrices 
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Then, it is easy to show that the concurrence becomes 

C12 = C u = C 23 = 2/. (14) 

Figure ((5]) shows the stationary entanglement vs 7 for differ- 
ent initial state. For initial states like \egg) 1 \geg), \gge) the 
entanglement does not depend on 7 . For |eee) initial state, 
by increasing the difference in the dissipation rates more en- 
tanglement will be induced in the system. On the contrary for 
\ege) the maximum value of entanglement is achieved when 
the dissipation rates into the two environments are the same. 
As it is expected the entanglement for initial states \eeg), \gee) 
are equal at 7 = 1/2. In these cases, the maximum value can 
be attained when there is one excitation in the initial state. 



IV. THE THREE QUBITS CASE WITH CLOSED 
BOUNDARY CONDITIONS 



We now study the chain of three sites with closed boundary 
condition. In this case the master equation (f2| reads 



dp 

~dt 



= 7 [2(^1 + ° 2 )p(vi+°t) 
-{(&+ + &+)(&- + &-) J}] 

+ P [2(ctJ +^)/5(<7 2 f +a£) 

-{(*3+*t)(&3+*r),p}] > 



(15) 



where the possibility for each qubit of having an asymmetric 
decay rate on the left and right environments is accounted for 
by the real factors 7, fj, and v (0 < 7, fi, v < 1). Clearly the 
symmetric situation is recovered when 7 = /i = v (= 1/2). 

Again arranging the density matrix (expressed in the com- 
putational basis {|e), |g}}® 3 ) as a vector v (e.g. writing 
Pi.j = v 8(i-i)+j)> me master equation (T5[ can be rewritten 
as a linear set of differential equations 

v(t) = Mv(i), (16) 

where M is a 64 x 64 matrix of constant coefficients given by 



M = 7(2ii (8> L-l - l\L\ <g> ig - I 8 ® l\L-l) 
+ fi(2L 2 ®L 2 - L ] 2 L 2 ®h-h® L\l 2 ) 
+ v(2L 3 <g)L 3 - LIL3 <g> I 8 - I s ® l\l 3 ), 



(17) 



where 

L\ = a~ ® I 2 ® I 2 + h 
L 2 = I 2 ®cr~®I 2 + I 2 

L 2 = I 2 ® I 2 ® cr _ + cr _ 

In this case it is possible to see that 
Mv = 0. 



er~ ® I 2l 
I 2 ®a~, 



(18) 



(19) 



admits a unique solution Vj = 5j 64> i- e - me on ly possible 
solution is \ggg), for all values of 7, fi, v. Hence no entangle- 
ment survive at stationary conditions. 



V. CONCLUSION 

In this work we have considered a spin-i chain dimerized 
by environments. By means of dissipative mechanism, each 
environment induces a chain link between contiguous qubits. 
Then we have studied the possibility of having long living en- 
tanglement without resorting to any other mechanism. In par- 
ticular for the case of three qubits chain with open boundary 
condition we have classified the amount of stationary entan- 
glement accordingly to some initial (separable) states. Here 
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we have also shown the appearance of the entanglement sud- 
den birth. On the contrary, for the case of three qubits chain 
with closed boundary condition we have proved the impossi- 
bility of stationary entanglement. This fact can be interpreted 
as entanglement frustration phenomenon ifiofl, induced in this 
context by the imposed periodic boundary conditions. 

The proposed scheme lends itself to to be extended to n > 
3 sites where one can evaluate how entanglement scales as 
function of the distance between the two qubits. Moreover, in 
the limit of large n, it could also results useful for studying 
possible connections between quantum phase transitions and 
reservoirs properties. 

Finally, the discussed method of generating entanglement 



seems economical and offers interesting perspectives for the 
generation of the so called graph states (useful for measure- 
ment based quantum compuatation) Intl . when one considers 
a network topologies more complicate than a simple chain. 
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